Point-Contact Conductances at the Quantum Hall Transition 
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On the basis of the Chalker-Coddington network model, a numerical and analytical study is made 
of the statistics of point-contact conductances for systems in the integer quantum Hall regime. In 
the Hall plateau region the point-contact conductances reflect strong localization of the electrons, 
while near the plateau transition they exhibit strong mesoscopic fluctuations. By mapping the 
network model on a supersymmetric vertex model with GL(2|2) symmetry, and postulating a two- 
point correlator in keeping with the rules of conformal field theory, we derive an explicit expression 
for the distribution of conductances at criticality. There is only one free parameter, the power law 
exponent of the typical conductance. Its value is computed numerically to be X t = 0.640±0.009. The 
predicted conductance distribution agrees well with the numerical data. For large distances between 
the two contacts, the distribution can be described by a multifractal spectrum solely determined by 
X t . Our results demonstrate that multifractality can show up in appropriate transport experiments. 

PACS numbers: 73.23.-b, 73.40.Hm, 61.43.Hv 
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I. INTRODUCTION 

Models of two-dimensional (2D) noninteracting elec- 
trons subject to disorder and a strong magnetic field, 
form what is called the (integer) quantum Hall universal- 
ity class. Their most prominent feature is the existence 
of a localization-delocalization (LD) transition, which un- 
derlies the plateau-to-plateau transition of the Hall con- 
ductance observed in the integer quantum Hall effect. 
Among the various members of the quantum Hall uni- 
versality class, the Chalker-Coddington network model 
jjj has been found |^|| to be a convenient representative, 
particularly for numerical purposes. A wavefunction in 
this model is a collection of complex amplitudes, one for 
each bond of a square lattice. The time evolution oper- 
ator acts on the wavefunctions by discrete steps, which 
are determined by unitary scattering matrices assigned 
to the vertices of the lattice. It has been established 
that, with varying left-right asymmetry of the scattering 
probability p € [0, 1], the stationary states undergo a LD 
transition with a critical exponent v w 2.35 for the 
localization length: £ ~ \p — p*\ ~ u : p* = 0.5. Moreover, 
it was shown that the critical states (for which the lo- 
calization length is much larger than the system size L) 
have multifractal properties || that are universal, and 
are characterized by an exponent ao ~ 2.28 describing 
the scaling of the typical value (i.e. the geometric mean) 
of the squared amplitude: exp (\n \ip\ 2 ) ~ L~ a °. The 
combination v(a.Q — d) is the critical exponent of the typ- 
ical local density of states (LDoS), which has been argued 
H to be an order parameter for the LD transition. 

The critical exponent v was extracted in a number of 
transport experiments In contrast, no such experi- 
ment has been carried out to measure ao or related mul- 



tifractal exponents. Multifractality is observable in local 
quantities such as local densities, or in the wave num- 
ber and frequency dependent dynamic structure factor 
(cf. HQ ). These have yet to be studied in transport 
measurements under mesoscopic conditions. In || it was 
pointed out that multifractality can show up in the size 
dependence of the conductance distribution close to the 
LD transition. Further suggestions for an experimental 
determination were made in ||, but it seems that these 
are still awaiting realization. Recently, it has been sug- 
gested that multifractality relates to the corrections to 
scaling and may be observable in the temperature de- 
pendence of the peak-hight of the conductance at the 
LD transition . In the present paper we demonstrate 
that another sensitive probe are point-contact conduc- 
tances, making multifractality directly accessible through 
a suitable transport measurement. By a point-contact 
conductance we mean a conductance between two small 
interior probes separated by a distance r. They show 
strong mesoscopic fluctuations at the LD transition, sim- 
ilar to those of the conductance (see e.g. []ll"|-|l4j| ) . By 
varying the distance r, point-contact conductances allow 
to study local details of mesoscopic fluctuations that are 
not captured by the (global) conductance. 

One of our motivations came from Ref. Jl5| , which 
discusses the conductance between two (or more) small 
interior contacts from a field theoretic perspective. In 
that work, it was pointed out that for a large enough 
system the point-contact conductance depends only on 
the distance r between the interior contacts, whereas the 
typical LDoS also involves the system size L. At criti- 
cality, the conductance in the infinite plane falls off alge- 
braically with r and, by the conformal hypothesis, this 
decay should be conformally related to the decay in other 
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geometries such as the cylinder. Unfortunately, a com- 
prehensive analytical theory of the critical point does not 
yet exist, and we cannot predict the critical exponents. 
As we are going to show, however, we are able to draw 
some strong conclusions just from the assumption of the 
existence of a conformal field theory for the critical point. 

In the present work, we use the Chalker-Coddington 
model to calculate point-contact conductances as well as 
quasi energies and the corresponding stationary states. 
In Sec. H we consider the dynamics of a closed network. 
To begin with familiar objects, we demonstrate that, un- 
der critical conditions, the local density of states has mul- 
tifractal correlation exponents that agree with those ob- 
tained from Hamiltonian models |lq] . These results will 
later be contrasted with the critical behavior of point- 
contact conductances in Sec. [V. We then review the 



definition of the point-contact conductances and com- 
pute them in the localized regime. We find them to be 
well described by a log-normal distribution determined 
by a sing le parameter, the typical localization length £ t 
(Sec. III). Our main theme, the investigation of point- 
contact conductances at p* = 1/2, is taken up in Sec. IV, 



which employs a combination of analytical and numeri- 
cal techniques, and is divided into five subsections. The 
q th moment of the conductance, (T q ), is expressed as a 
two-point correlator of a supersymmetric vertex model 
equivalent to the network model |l5|,[l7],[l8| . We exploit 
the global GL(2|2) symmetry of that model and, based 
on the general principles of conformal field theory, pro- 
pose a closed analytical expression for (T q ). From that 
we extract predictions for the typical conductance, the 
log-variance and the entire distribution function. These 
predictions leave one parameter undetermined, which can 
be identified with the power law exponent X t of the typ- 
ical conductance. We calculate X t numerically and find 
X t = 0.640 ± 0.009. The distribution function with this 
value is shown to agree with the numerical data. In the 
limit of large separation between the contacts, which is 
difficult to reach numerically, the conductance distribu- 
tion can be reduced to a spectrum of multifractal expo- 
nents. We speculate on a possible connection between the 
multifractal spectra of the point-contact conductances 
and of the local density of states. 

Finally, by matching to results that are available for 
the quasi-lD limit, we argue that, if the network model 
(or a suitable continuum limit thereof) were a fixed point 
of the renormalization group, then the power law ex- 
ponent of the typical conductance would have to be 
X t = 2/tt ps 0.637, which is in remarkably close agree- 
ment with our result from numerics. Because standard 
conformal field theories of the Wess-Zumino-Witten or 
coset type predict critical exponents to be rational num- 
bers (if the fixed point is isolated) such a value, if correct, 
would imply an unconventional fixed point theory. 



II. MULTIFRACTALITY OF THE LOCAL 
DENSITY OF STATES AT CRITIC ALITY 

To describe the critical behavior of 2D electrons in a 
strong magnetic field and a smooth random potential, 
Chalker and Coddington formulated a network model 
which is very simple and yet captures the essential fea- 
tures. The model is composed of a set of elementary 
"scatterers" placed on the vertices of a square lattice. 
(In a microscopic picture, these correspond to the saddle 
points of the random potential p9[.) 




FIG. 1. Graphical representation of the Chalker-Codding- 
ton network model. Wave amplitudes propagating on the 
links are scattered to the left (right) with probability p (resp. 
1 — p). The corresponding unitary scattering matrices S are 
situated at the nodes of the network. 

Unidirectional channels link the scatterers to each 
other as shown in Fig. 1. Each elementary scatterer is 
represented by a unitary 2x2 scattering matrix S trans- 
forming 2 incoming amplitudes into 2 outgoing ampli- 
tudes. To simulate the effect of a disordered or irregular 
array of scatterers, a (kinetic) phase factor is attached to 
every channel amplitude. These phases are taken to be 
independent random variables distributed uniformly on 
the interval [0, 2ii\. An incoming amplitude on a given 
link can only be scattered to the left or right. Let us de- 
note the probability for scattering to the left by p. The 
scattering probability to the right is then 1 — p. The pa- 
rameter p for each scatterer can be taken to be fixed or 
drawn at random from a certain distribution. In either 
case, the states in a system of finite size L turn out to 
be delocalized jl],|| when the mean of p lies in a small 
interval around p = 0.5. The width of this interval, Ap, 
shrinks to zero in the thermodynamic limit, Ap ~ L^ 1 ^ ', 
where v « 2.35 is the critical exponent of the localization 
length £. 

A network model wavefunction is a set of complex 
amplitudes {ip(l)} where / = l,...,Ni runs over the 
links of the network, and the normalization is fixed by 
W)| 2 = 1- Wavefunctions are propagated forward 
in time by discrete steps [p0|H , 
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(i) 



i'=i 



the propagator for one unit of time being a sparse unitary 
Ni x Ni matrix U which is uniquely determined by the 
network of scattering matrices |2lj . 

Stationary states of the (isolated) network are solu- 
tions of the equation PJ2^| 
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The eigenphases <p n will be referred to as quasi ener- 
gies. (This terminology is motivated by the analogy with 
the eigenphases of the Floquet operator of a periodically 
driven quantum system.) The stationary states of the 
network model are critical (which is to say their local- 
ization lengths £ exceed the system size L) when the pa- 
rameter p is close enough to p* = 0.5. The critical states 
are multifractals Criticality is visible through power 
law scaling of the moments §|: (\ip(l)\ 2q ) L ~ L~ d - T ^ 
where r(q) is a nonlinear function of q, and d = 2. (We 
will continue to write d instead of 2 to keep the depen- 
dence on the number of dimensions explicit.) Within 
our numerical accuracy, the r(q) functions for different 
critical states coincide. The corresponding distribution 
function is specified by a single-humped positive function 
/(a), called the multifractal spectrum of the wavefunc- 
tion @, 

prob(F = |?A(/)| 2 ; L)dP ~ L~ d+f ^da , 

where InP = — ahiL, and f(a) is related to r(q) by a 
Legendre transformation f(a(q)) = qa(q) — r(q), a{q) = 
dr{q)/dq. In the vicinity of its maximum, /(a) is well 
approximated by a parabola p3] : 



/(«) 



(a - Q ) 2 
4(a - d) ' 



which is seen to be determined by a single number olq. 
From jl6|,[3| we know this number to be «o ~ 2.28. 

The local density of states (LDoS) is defined as 
p(4>J) — ^ n 5{4> ~ ^n)\fa-{l)\ 2 where fa denotes a sta- 
tionary state with quasi energy <fr n . To wash out the 5- 
peak structure of this function in a closed finite system, 
we smoothen it over a scale of one mean level spacing, 5. 
The LDoS then becomes 

p(cj>,l) = 5- 1 \i;(<p,l)\ 2 , 

where \ij>(<p, l)\ 2 means the square of the wavefunction 
amplitude, microcanonically averaged over the quasi en- 
ergy interval [4> — 5/2, <j> + 5/2]. Given 5" 1 ~ L d and the 
multifractal scaling law for the critical states, the LDoS 
must scale as (p(4>,l) q ) L ~ L~ Ap ^ with 

A p (q) = (l-q)d + r(q) , (2) 

and the typical value as p t = exp(lnp(l)} L ~ L d ~ a ". 



The average LDoS is well known to be noncritical and 
nonvanishing at the LD transition. In contrast, the typ- 
ical LDoS does show critical behavior. For one thing, it 
is zero in the region of localized states. For another, sup- 
pose the system under consideration had a finite band 
of metallic states (which it does not), as is the case for 
a mobility edge in three dimensions. The typical LDoS 
would then be nonzero in that metallic band and, for 
values of L much larger than the correlation length £ c , 
would vanish with exponent /3 t = ^(«o — d) on approach- 
ing the critical point. Such behavior is reminiscent of an 
order parameter, which distinguishes between two phases 
joined by a second order phase transition. The LDoS has 
in fact been proposed as an order parameter field for 
the general class of LD transitions Q . Although the 2D 
quantum Hall class has no extended metallic phase, the 
exponent ao — 2m 0.28 still controls the finite size scaling 
of the LDoS close to the critical point. 

It is then natural to inquire into the nature of the cor- 
relation functions of the LDoS at criticality. This was 
done in |p5 16 (see also (26)). The critical correlations 
turned out to be of the form 



p [ $ (u},r = \r\,L) = (p(e, l) q p(e + u,l + t)*) 1 
~ {r/L^Y^L-^ , 



(3) 



where the length L u is defined as the linear size of a sys- 
tem with level spacing u: — (uj/6)~ 1 / d L ||. This 
length scale provides a natural cutoff for the correlations 
between two critical states ||2~7| , ^ , |l6| ,[7| , whence @ ap- 
plies to the regime r -C L w < L. For L u > L the depen- 
dence on L u saturates and L w has to be replaced by L in 
([jj). The critical exponent z{q) is given by the important 
scaling relation p9|,p5[| 



z(q) = d + 2r{q)-T(2q). 
A second scaling relation, 

%) = 2(1 -q)d+2 T (q) = 2A p (q) 



(4) 



(5) 



follows from (g) and (|J) by letting r and ui go to zero 
and then matching with {p 2q )h ~ L~ Ap ( 2q \ We empha- 
size that the multifractal correlations of the LDoS have 
a characteristic dependence on system size Q, which is 
encapsulated by the family of exponents z{q). 

We now report on a numerical test of the scaling rela- 
tion for z(q), using the Chalker-Coddington model. Let 
us start by verifying that the critical correlations are in- 
deed cut off at the scale L u . 

Fig. 2 plots —ln(p^\uj,r, L) as a function of lnr for 
q = 0.5, L = 150, and two different values of ui. For 
u> = 0.000541 (corresponding to w 76) the linear de- 
pendence on In r is seen to be limited to a much smaller 
region than for lo = 0.00026 (L u « 110). The rest of a 
our data show similar behavior, although the fluctuations 
due to finiteness of the data set are quite large. It is clear 
though from these data that does set the typical scale 
for the cutoff of critical correlations. 
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FIG. 2. Double logarithmic plot of the correlation function 
P2° J ' (w, r,L — 150) versus distance r for two different energy 
separations u>, which correspond to cutoff scales L^. 



To study the scaling exponent for the r-dependence of 
the correlator, we first took L u — oo, i.e. zero energy 
separation, and compared the calculated exponent z(q) 
with the value predicted by the scaling relation (|4|) and 
the r(q) spectrum previously obtained. As shown in the 
upper part of Fig. 3, the result for z(q) agrees with the 
prediction within the statistical errors. Next, we took 
the quasi energies to be separated by a finite amount 
corresponding to L u ss 110, and calculated z(g) for a few 
selected values of q. Within the error bars, agreement 
with the previous values was obtained. 

Finally, we investigated the scaling with respect to L u 
which, according to (Q), should lead to the same z(q). 
The data set consisted of 500 eigenfunctions for a par- 
ticular disorder realization and a range of quasi energies 
permitting small values of L w to be reached. The system 
size was L = 50 and we performed a spatial average. The 
choice of a relatively small system size was necessitated 
by the fact that the fluctuations of the eigenfunction cor- 
relations grow in strength as lo is increased. To reach 
statistical convergence, a very large number of pairings 
p(4>, l)p{4> + oj, I + r) for fixed values of r and u) must be 
accumulated. This, given present computer capacity, is 
possible only for a small enough system. We set r = 2 
and found behavior sufficiently linear in In L u in a regime 
between = 3 and L u = 12. The results for z(q) are 
shown at the bottom of Fig. 3. In view of the systematic 
difficulties, we conclude that our data are consistent with 
the scaling relation. 

In summary, our numerical results for the multifractal 
exponents of the local density of states clearly identify the 
Chalker-Coddington model as belonging to the quantum 
Hall universality class. 




FIG. 3. Scaling exponents z(q) obtained from scaling with 
respect to distance r (top) and energy separation ui (bottom). 
The data points denoted as "numerics" correspond to the lin- 
ear regression of In versus In r (top) and In (bottom) 
while the data denoted as "scaling" follow from the scaling 
relation (jij) and the known r(q) spectrum. 



III. POINT-CONTACT CONDUCTANCES IN 
THE LOCALIZED REGIME 

We start by reviewing what is meant by a "point- 
contact conductance" in the Chalker-Coddington model 
p^| , p7t . Select two of the links of the network for "con- 
tacts" and cut them in halves (Fig. 4). Then do the fol- 
lowing. Inject a total of one unit of probability flux into 
the incoming contact links, and apply the time evolution 
operator U once. Then inject another unit of flux and 
apply U again. Keep iterating this process, always feed- 
ing in the same unit of probability flux so as to maintain 
a constant current flow into the network. The outgo- 
ing ends of the broken links serve as drains, so flux will 
eventually start exiting through them. After sufficiently 
many iterations of the procedure, the network will have 
settled down to a stationary state. The stationary wave 
amplitudes at the severed links square to transmission 
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and reflection probabilities, which translate into a (point- 
contact) conductance by the Landauer-Biittiker formula. 



contact 1 




contact 2 



FIG. 4. A section of the network model with two interior 
point contacts. 

In the present section we describe how to obtain these 
amplitudes by solving a system of linear equations in- 
stead of laboring through many iterations of the above 
dynamical procedure. Afterwards, we calculate the dis- 
tribution of point-contact conductances in the localized 
regime and show that it is well approximated by a one- 
parameter family of functions, depending only on the 
value of the typical localization length £ t ■ 

The injection of current into the contact links c, d is 
simulated by modifying the time step in the following 
way: 

KW = U(\iIh) + a\c) + b\d)) , 

where a and b (subject to \a\ 2 + \b\ 2 = 1) are the ampli- 
tudes of the current fed into the links c and c', and |c), \d) 
denote basis states with unit amplitude at c, c', respec- 
tively, and zero elsewhere. To implement the draining 
action of the outgoing ends at c and c', we define pro- 
jection operators Pc by Pc\"^) — ip(C)\C) for C — c,c'. 
Using these, we can write the complete dynamics as 

hk=o) = o, 

|^t + i) = U((l - P c )(l - P c ,)|Vt) + a\c) + b\c')) . 

The stationary current carrying state is formally ob- 
tained by taking the limit \ipoo) — liui{— «x> \ipt)- Alterna- 
tively, we set Up = {7(1 — P c — P c <) and use stationarity 
to deduce for \ipoo) the linear equation 

(l-CiOWU = U(a\c)+b\d)), (6) 

which can be solved by an inversion routine. 

What we want are the amplitudes of the stationary 
state at the links c and c', which are the components 
V^oo(c) and V'oo(c') of the vector |V>oo)- According to the 
Landauer-Biittiker formula, the conductance g between 
two point contacts c and d is given by the transmission 
probability T = \t c > c \ 2 as g — (e 2 /h)T. The scattering 



problem for this is defined by feeding one unit of current 
into the link c and zero current into d . The transmission 
amplitude is the amplitude of \ipoo) at the exit d. Hence 
we set a = 1 and 6 = 0, and compute the conductance 
from 

r = |(c'|^oo>| 2 . (7) 

Note that such a point-contact conductance is bounded 
from above by unity: T < 1. 

Our numerical calculations of point-contact conduc- 
tances were done for systems of size L = 60, 80, and 100, 
with the distance between contacts varying from r = 1 
to r = L/2. For every distance we calculated between 
200 and 3000 conductances (depending on system size). 
In the following, we present results for the distribution of 
the conductance for localized states ((<L), correspond- 
ing to the plateau regime of the quantum Hall effect. 
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FIG. 5. Normalized histogram of the logarithm of the 
conductance T in the localized regime, for a network of linear 
size L = 60 and a fixed distance r = 7 between the contacts. 

Fig. 5 shows a normalized histogram accumulated from 
3000 data points for the logarithm of the conductance 
g = T (in atomic units), at fixed values of the distance 
r = 7 and the system size L = 60. The histogram 
clearly demonstrates the approximate log-normal char- 
acter of the distribution. Such behavior is expected from 
the standard picture of localization. More precisely, the 
picture says that the inverse of localization length £ has 
a normal distribution (see pj for a review), and the con- 
ductance is given by g = go exp(— 2r/£) where the factor 
of 2 is due to the convention of associating the localiza- 
tion length with the modulus of the wavefunction. The 
inverse of the average of l/£ is called the typical value 
£t of the localization length: = (C -1 )- As a further 
check on this picture, we calculated (In T) as a function of 
r and found linear behavior, as shown in Fig. 6. The av- 
erage was performed as a combination of spatial average 
and disorder average accumulating 3000 data. 

A more interesting question concerns the relation be- 
tween the two parameters (log-mean (InT) and log- 
variance ((<51nT) 2 )) of the log-normal distribution. Ac- 
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FIG. 6. Plot of the average of InT versus the distance r 
between the contacts, for different system sizes L and proba- 
bilities p = (1 + e~ E y 1 . 



cording to the one-parameter scaling hypothesis first for- 
mulated in p0| , one expects the two parameters to be 
dependent on each other. Such a dependence was indeed 
observed in other LD transitions (see e.g. and it 

takes the form of a linear relation 



((c51nT) 2 ) = -A(\nT)+B. 



(8) 



where A is a number of order unity (it equals 2 in quasi- 
1D systems fl32"[), and B is some offset due to the presence 
of the factor g Q . Our data confirm this picture. The 
value for A we find is A = 1.00 ±0.05 (Fig. 7). Thus, the 
typical localization length is the only relevant parameter 
for the conductance distribution in the localized regime. 
This distribution is very broad and is well described by 
a log-normal form. 
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FIG. 7. Linear relation between the mean and the variance 
of the conductance distribution in the localized regime. 



IV. POINT-CONTACT CONDUCTANCES AT 
CRITIC ALITY 

We now embark on an investigation of point-contact 
conductances of the critical network model at p* = 1/2, 
by employing a combination of ana lytical and numeri- 
cal techniques. In subsection IV A, the network model 



is mapped on a supersymmetric vertex model (for tech- 
nically related work see also |1S[] ) and the q th moment 
of the transmission coefficient T is expressed as a two- 
point correlator (Eqs. 
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i 12|)). For this we follow the 
method of fTij] , which is based on the so-called color- 
flavor transformati on an d is included here for complete- 
ness. In subsection IV B we take advantage of the global 
symmetries of the vertex model and, based on the general 
principles of conformal field theory, propose a closed ana- 
lytical expression for the <; th moment of the point-contact 
conductance at criticality (Eq. (|l7|)). This expression is 
analytically continued to q — 0, to extract predictions for 
the typical conductance and the log- varian ce, w hich are 
compared to numerical data in subsection IV C . In the 
limit of large separation between the contacts, which is 
difficult to reach numerically, a multifractal description 
of conduct ances becomes appr opria te, as is discussed in 
subsection [VP . In subsection IV E] we use conformal in- 
variance to predict an exact expression for point-contact 
conductances in the quasi one-dimensional limit of the 
network model. Finally, in subsection IV F we recon- 



struct from the moments the entire distribution function 
of point-contact conductances. 



A. Mapping on a supersymmetric vertex model 



We start by recalling the Landauer-Biittiker formula 
for the dimensionless conductance: g = T = \t c > c \ 2 , where 

td e = (c'Ki-Upy'uic) . 

For the following it is convenient to change our conven- 
tions slightly and interpret each of the severed contact 
links c and d as a pair of disjoint links (Fig. 4), defining 
basis states |oi) and |o 2 ) respectively. We then 
impose the boundary conditions 



U\ 0l ) = = U\o 2 ) , 

u^\h) = o = trt|i 2 ) 



(9) 
(10) 



The boundary conditions on the outgoing links have the 
same effect as the projectors P c and P c r, which enter in 
the definition Up = U(l — P c — P c ')- With these bound- 
ary conditions in force, we can write the expression for 
the transmission amplitude as t c > c = (02KI — t/) 1 1 « 1 } - 
Assuming the two point contacts to be separated by more 
than one lattice unit, so that (o2\U\ii) — 0, we can write 
it in the even simpler form 



tr 



(021(1-11)-%) 
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Note that the escape of flux, here modeled by the bound- 
ary conditions on the outgoing links, causes a unitar- 
ity deficit and thus ensures positivity of the operator 

1 - Wu. 

The next step is to express t c i c as a Gaussian su- 
perintegral. For this purpose we introduce a doublet 
^ = (^Bj^f) for every link. (The incoming and out- 
going ones are included.) The quantities 'J'b and \I/f are 
bosonic and fermionic (i.e. commuting and anticommut- 
ing) complex integration variables. Defining a Gaussian 
"statistical" average by 



exp —5+ 



with 



s+= E f -(O(^-wM0, 



er— B,F L, 



we have 



tc'c = <*b(o 2 )*b(«i)> • 

As usual the integration measure for the field "J" is the fiat 
one, normalized so that (f) = 1. The Gaussian integral 
over <3> converges because the modulus of the operator U 
with boundary condition (^) is less than unity. 

Our interest here is not only in the average conduc- 
tance but in all moments of the conductance. In fact, we 
ultimately will reconstruct the entire distribution func- 
tion. To calculate the q moment we need an expression 
for the q th power of t c > c . As is easily verified from Wick's 
theorem, this is given by 

4 = ^B( 02 )%( ? ir) . 

Note that in a many-channel situation (such as a network 
model with more than one channel per link or a conduc- 
tance which is not point contact) where t c > c consists of 
more than one amplitude, going from q = 1 to arbitrary 
q requires enlarging ^ to a superfield with more than one 
component. Here we are in the fortunate situation that 
a single component suffices to generate all the moments 
II- 

We call <3/ a "retarded" field and denote it from now 
on by V E'+. The complex conjugate t c / c is expressed by 
a similar construction using an "advanced" field, 
By combining the Gaussian integrals over retarded and 
advanced fields we get 



jig 



1 



(* +b (o 2 )^ +b (*i) 9 *-b(h) 9 *-b(o 2 )' ? ) 



for the g th moment of the transmission probability. The 
Gaussian statistical average here is taken with respect to 
exp — S where S = S+ + S— and 

S-= Y,*-*(l)(Sii'-Uin)*-„(l') ■ 

<T=B,F 1,1' 



Our next goal is to average over the disorder of the 
network model. For that we put 

U vl = U x {l',l) e^W 

where U\ is the deterministic part of U describing the 
scattering at the nodes, and ip(l) are random phases uni- 
formly distributed on the interval [0, 2iz\. On extract- 
ing from exp— S the (^-dependent parts and temporarily 
omitting the integration over the superfield 'J, we are 
faced with the integral 

J l 77 <j;l',l 

+ ^^{l)^ Ml) U\{l,l')^-a{l') 



It is seen that the random phases e ±tv couple, roughly 
speaking, to the bilinears \J r ± (T \f r ± ( T. Plain integration 
over the variables ip(l) now produces a product of Bessel 
functions, an expression which is not a good starting 
point for further analysis. Fortunately, there exists some- 
thing better we can do. In the supcrsymmctric treat- 
ment p4| of disordered Hamiltonian systems, one makes 
a Hubbard-Stratonovitch transformation replacing the 
disorder average by an integral over a supermatrix- valued 
field Q. It turns out that a similar replacement, called 
the "color-flavor transformation" |55|, can be made in 
the present context. The role of Q is taken by two sets 
of complex fields I 1— > Z aa >(l) and I 1— > Z^^il), which 
assemble into supermatrices, 

Z = ( Zbb Zbf I Z = ( '? BB Z ~ B¥ \ 

y Zfb Zpp J ' y Zfb Zpp J 

and couple to 4 r + (T \f r _ cr / and ^ _ (T / 1 S' +cr respectively. The 
final outcome of the transformation will be another ex- 
pression for given at the end of the next paragraph, 
where instead of integrating over the random phases <p(l) 
we integrate over the supermatrix fields Z, Z. 

For completeness, let us now briefly summarize the 
mathematical structures surrounding the supermatrices 
Z and Z. The details can be found in |3(|. Elements 

g = of the Lie supergroup Gc = GL(2|2) are 

taken to act on Z by the transformations 

Z h-> g ■ Z = (AZ + B){CZ + D)- 1 , 
Z m. g ■ Z = (C+ DZ)(A + BZ)- 1 . 

Because this group action is transitive and Z = Z = 
is fixed by an He = GL(1|1) x GL(1|1) subgroup with 

elements h = ( ^ D^j ' ^ 6 com P^ ex supermatrices 

Z and Z parameterize the complex coset superspace 
Gc/ He- Let D(Z,Z) be a Gc invariant superintegra- 
tion measure on this coset space and put D/j(Z, Z) = 
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D(Z, Z) SDet(l — ZZ). We take the integration domain 
for the bosonic variables to be the Riemannian subman- 
ifold M r defined by the conditions 

|2 



Zpp — — z 



FF 



■2bb — +Z\ 



BB 



J BB | 



< 1 



The normalization of D (Z, Z) is fixed by requiring that 
f Dfi(Z, Z) = 1. Given these definitions, it was shown in 
] 36| that the above expression for J- is equivalent to 



T= I T>n(Z,Z) exp^ (*_ CT ,(0-W0*+ CT (0 

+(§ +(T c/ 1 )(/)z^(0(t/ 1 t *- CT 0(0 



where V/j,(Z,Z) = H l Dfx(Z(l), Z{1)). It is understood 
that the sum over / excludes all links leaving the net- 
work. Indeed, by the boundary condition (||) the phase 
ip(l) on an outgoing link I = o\ or I = o 2 never appears 
in the formalism and there is no need to introduce any 
supermatrices Z, Z there. Alternatively, we may extend 
the sum to run over all links and compensate by putting 

Z(l) = Z(l) = for all outgoing links. 

We mention in passing that the equivalence between 
the two expressions for T extends |17]| to the more gen- 
eral case where the random phase factors e lip ^ € U(l) 
on links are replaced by random U(N) elements. The 
transformation from these U(7V) "gauge" degrees of free- 
dom to composite "gauge singlets" Z(l),Z(l) is in strik- 
ing analogy with strongly coupled U(iV) lattice quantum 
chromodynamics, where integration over the gauge fields, 
which carry U(7V) color, produces an effective description 
in terms of meson fields carrying flavor. Therefore, the 
passage from the first form of Z to the second one is 
referred to as the "color-flavor transformation" |35| . 

An attractive feature of the color-flavor transformation 
is that it preserves the Gaussian dependence of the inte- 
grand on the fields "J, ^ . If we switch to schematic nota- 
tion and suppress the link and super indices, the Gaus- 
sian statistical weight is the exponential of the quadratic 
form 



Z 



-UyZU\ 
1 



In the absence of source terms or other perturbations, 
integration over ^> simply produces the inverse of a 
superdeterminant , 



SDet~ 



1 

-Z 



-u y zu\ \ _ 



= SDet _1 (l - ZU X ZU\) . 



When the source terms are taken into account, addi- 
tional factors arise. Since we are calculating a point- 
contact conductance, these are localized at two points. 
By applying Wick's theorem, we get from the combina- 
tion 4 , _b(«i) 9 ^ , +b(«i) 9 /9! at the first contact an extra 
factor 



Z(l - U X ZU\Z) 



bb 



(ii,h) q = Z^ B (h) 



The simplification to the right-hand side of this equation 
occurs because Xj\ — according to the boundary 
conditions (10). Unfortunately, a similar simplification 
does not take place at the other contact. (There is a 
basic asymmetry between "in" and "out" , which stems 
from the fact that we must decide on some ordering of the 
factors e llp and U\ in the evolution operator for one time 
step.) However, we can mend the situation with a little 
trick. We "prolong the exit" by inserting an additional 
node on the out-links as shown in Fig. 8. 




FIG. 8. Prolongation of the exit by adding a dummy 
vertex at the drain contact, so as to simplify the field theoretic 
representation (see text) of the point-contact conductance. 

This means that flux arriving at the link 02 does not 
exit immediately, but first gets transferred to o 2 and then 
exits from there during the next time step. With this 
modification, the extra factor arising from the product 
y + B(o' 2 ) q ^- B (o' 2 ) q /q\ simplifies to 



u x zu\{\ - zu x zul)- x (o' 2 ,o' 2 y 

' BB 



since U\{o 2 ,o 2 ) = U\(o 2 ,o' 2 ) = 1 and Z(o' 2 ) = 0. 

In summary, we have expressed the q th moment of the 
transmission coefficient T as a two-point correlator, 

<T«) = (z^ B ( 02 )^ B (n)) , (11) 

of a supersymmetric lattice field theory, 

(•) = J T>fjt(Z,Z) • SDet _1 (l - ZUyZUl) . (12) 

The superdeterminant runs over superspace and the 
Hilbert space of the network model. 

As it stands, the formula for (T q ) does not display 
clearly the internal symmetries of the theory. To make 
these more explicit, recall the definition of the (single Z) 
integration measure, 

Dfi(Z, Z) = D(Z, Z) SDet(l - ZZ) . 

We now factor the superdeterminant as two square roots, 

SDet(l - ZZ) = SDct 1/2 (l - ZZ) x SDet 1/2 (l - ZZ) , 



^bb(°2) , 
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and associate each of these with one of the two nodes a 
link begins or ends on. This procedure works perfectly 
for the internal links: every square root factor is assigned 
to exactly one node. For the external (or contact) links 
the situation is different. These are connected to just 
a single node of the network, and therefore only one of 
the two square roots is used up. It is natural to assign 
the remaining factor to the operator whose correlation 
function we are evaluating. In this way, the expression 
for (T q ) reorganizes to 

(T q ) = (v q (o 2 )v* (i 2 ) x vo(oi)vl q (ii)) (13) 

where 

v q = Z^ B SDet 1/2 (l - ZZ) , 
v - q = ^bb SDct 1/2 (l-ZZ) , 

and 

(•) = J T>(Z,Z) • exp-S , 

V(Z,Z)=l[D(Z(l),Z(l)) . 
i 

The statistical weight exp —S is given by the superde- 
terminant SDet _1 (l — ZU 1 ZUl) multiplied by one (two) 
factors of SDet(l — Z(l)Z(l)) for every external (internal) 
link. As was shown in pM, exp— S has the structure of 
a vertex model weight, i.e. it separates into a product of 
factors, one for each node. Hence the lattice field theory 
with statistical weight exp — <S is called a vertex model. 
The factor for a given node, the so-called 1Z matrix, is 
invariant under the global action of GL(2|2) or, more 
precisely speaking, can be regarded jlTj as the matrix 
element of a GL(2|2) invariant operator between coher- 
ent states parametrized by the field variables Z, Z on the 
links which emanate from that node. 



B. Symmetries and conformal hypothesis 

The special significance of the functions v q and v*_ q 
is that they lie [|l5| in irreducible representation spaces, 
denoted by V and V* , of the symmetry group GL(2|2). 
The action of GL(2|2) on these spaces can be shown to 
be unitary for a SU(1, 1) x SU(2) subgroup, and V and 
V* belong to the discrete series of SU(1,1). The repre- 
sentation space V {V*) is of lowest weight (resp. highest 
weight) type. 

To make further progress, it is imperative that we 
exploit the global GL(2|2) symmetry of exp— S. The 
links i\ and 0\ lie close in space, and so do i 2 and o 2l 
whereas these two sets are in general far apart from each 
other. Therefore, our next step is to fuse v q (o 2 ) e V 
with VQ(i 2 ) £ V* and to make a decomposition into irre- 
ducible representation spaces of GL(2|2). The operators 
at the other contact, Uo(oi) and v*_ q (i\), are processed in 



the same way. Thus we need to know how to reduce the 
tensor product V ®V* . According to Sec. 5.2 of fL5|] , this 
reduction involves a single continuous series of GL(2|2), 
which is closely related to the principal continuous series 
of unitary representations of SU(1, 1). This series is la- 
beled by a real parameter A £ [0, oo). Denoting the basis 
state with weight m by ip\ m , we have 

Vo = J (Vq, V*0\\q) <p Xq [J,(X)dX , (14) 

V V*_ q = J (VO, V*-q\X-q) <p x - q fJ,(X)dX , (15) 

where (Vm, V*m'\Xm + m!) is a Clebsch-Gordan coef- 
ficient, and [i(\)d\ is the Plancherel measure for the 
continuous series labeled by A. Explicit expressions for 
these will be given below. We adopt the convention 
Tp\q~ = ipx—q. Then an immediate statement is that the 
Clebsch-Gordan coefficient behaves under complex con- 
jugation as 

(Vq,V*0\Xq) = (V0,V*-q\X-q) , 

as follows from v q Vg = v v*L q . 

All steps so far have been exact and rigorously jus- 
tified. Now we have to make an assumption which is 
crucial, namely that the vertex model at p* = 1/2 flows 
under renormalization to a conformal fixed point theory. 
Sadly, although a substantial effort has been expended on 
identifying that fixed point, we still do not understand its 
precise nature. Nevertheless, as we shall see, we can draw 
a number of strong conclusions just from the assumption 
of its existence. 

By the principles of conformal field theory, the two- 
point function of (p\ decays algebraically, 

(^(rVv- 9 (r')> = S -^jT^ 1 \r r']^ +■■■ (16) 

where A\ is the scaling dimension of the most relevant 
conformal field contained in the expansion of tp\q. (By 
GL(2|2) invariance, this dimension is independent of the 
weight q.) Although A a cannot be predicted without 
knowing the stress-energy tensor of the conformal field 
theory, we will make an informed guess later on. What 
we can say right away is that must be an even func- 
tion of A. (This is due to the invariance of A\ under 
a Weyl group action on the roots of the Lie superalge- 
bra of GL(2|2), which takes A into —A.) The appearance 
of the weight function /Lt(A) in the denominator, and of 
the Dirac 5-function (instead of the usual Kronecker 5- 
symbol), are forced by the fact that we are dealing with 
a continuous series. 

On inserting the decompositions (|l4|) and ( |l5| ) into the 
formula ( p"3| ) for (T q ), and using the expression ( fl6| ) for 
the two-point correlator, we obtain 

(T*) = / I {Vq, V*0\Xq) |V 2A >(A)dA , (17) 
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where r is the distance between the two point contacts, 
measured in the units that are prescribed by the choice 
of normalization made in ([l6]) . 

The computation of the Clebsch-Gordan coefficient 
and the Plancherel measure entering in the expression 
for (T q ) is nontrivial. Major complications arise from the 
fact that the modules V and V* are infinite-dimensional, 
and that the representation spaces appearing in the de- 
composition of the tensor product V ® V* have neither 
a highest nor a lowest weight vector. Consequently, the 
computation cannot be done solely by algebraic means 
and a certain amount of analysis must be invested. We 
have relegated this lengthy calculation to the appendix, 
where it is shown that 



H(X)d\ = -tanh 



dX 



(Vq, V*0\Xq) 



r(g 



f) r (<r 



1 ' 



r(# 



(18) 
(19) 



As it stands, the result for (T q ) has been derived for all 
positive integers q. However, since the Clebsch-Gordan 
coefficient is a meromorphic function of q, it is clear that 
we can analytically continue the result to all q. (This an- 
alytic continuation is unique, as follows from the bound 
\(T q )\ < 1 for Rc(q) > and Carlson's theorem stated 
in paragraph 5.81 of ]37].) Notice, now, that all of the 
g-dependence of (T q ) resides in the Clebsch-Gordan coef- 
ficient, while the dependence on r is encoded in the factor 
r -2A A would therefore seem that all moments of T de- 
cay asymptotically with the same power, (T q ) ~ r -2A ^ 
the exponent being given by Ao = Min^gR+AA. While 
this is correct for q > 1/2, it fails to be true for q < 1/2. 
The reason is that, when q is lowered past the value 1/2 
from above, the two poles at A = ±i(2q— 1) of the gamma 
functions T(q— ii-y) cross the integration axis at A = 0. 
Therefore, to do the analytic continuation to q < 1/2 cor- 
rectly, we must pick up the contribution from these poles. 
It is straightforward to calculate the residues at the poles, 
which results in 



(T q ) = T(qY 



|r(g- 



-2A, 



fj,(X)dX 



+ 27rcot(q7r)M r- 2A *(^> (\q\ < 1/2) . (20) 

Here we have used A>, = A_a to combine terms. 

By letting q go to zero in (^) and using the fact 
that normalization of the distribution function for T 
implies (T 9 )| _ Q = 1, we can deduce a constraint on 
A\. Indeed, from fli9| ) the square of the Clebsch- 
Gordan coefficient vanishes uniformly in A as q — > 0, 
and since lim g _>o 27rcot(g7r)r(2q)/r(q) 2 = 1, we have 



(T°) = 



Hence we conclude 
A_ 4 = A +i = 



Thus, normalization constrains A\ to be of the form 
A X = (X 2 + 1)F(X 2 ) . 



It will be convenient to set ^(A 2 ) = /(A 2 + 1). 



C. Typical conductance and log- variance 

Our next goal is to obtain information about the un- 
known function /. Because we have no analytical control 
on the fixed point theory, we are forced to resort to nu- 
merical means. To prepare the numerical calculation, we 
first show that /(0) is determined in a very simple way by 
how the typical conductance exp(lnT) varies with r. To 
that end we differentiate both sides of ( pp| ) with respect 
to q at q — and use the identity 



— (T q )\ 
dq K > lq =° 



(InT) 



The integral on the right-hand side of ([20|) tends to a 
finite value while T(q)~ 2 behaves as q 2 in the limit q — » 0. 
Also, the Taylor expansion of 27rcot(<77r)r(2<7)/r(<7) 2 = 
l + 0(q 2 ) contains no term linear in q, as is easily verified 
from standard properties of the gamma function. Hence 



(In T) = _ r - 2A «C2 9 -i) 



dq 

or, after exponentiation, 



9=0 



-8/(0) lnr 



exp(lnT) = 



, A t = 8/(0) 



(21) 



Thus, exp(lnT) decays as a pure power. Note that this 
feature is unique to the typical conductance. For a gen- 
eral value of g, formula (EOT) shows that the behavior of 
(T q ) as a function of r is governed by a whole continuum 
of exponents. 

We have calculated conductances at criticality (p 
>)' = 0.5) for systems of size L = 40, 60, and 100, and 
for distances r varying from r = 1 to r = L/2. For every 
distance, 2200 realizations for L = 40, 1200 realizations 
for L = 60, and between 200 and 300 realizations for 
L = 100, were generated. 

Fig. 9 shows 1300 realizations of the conductance for 
L = 60 and r = 10. The distribution of conductances is 
seen to be very broad. In Fig. 10 the mean logarithm of 
the conductance, (InT), is displayed as a function of lnr 
for three different system sizes L = 40, 60, 100. 

The typical value indeed scales as a power with r, 



cxp(lnT) = (r/r a )- x * , 

A t = 0.640 ± 0.009 . 



r w 0.42 , 



(22) 



and has no significant dependence on the system size. 
The result for X t is based on data for L — 100 and dis- 
tances r < 40, on data for L = 60 and distances r < 25, 
and on data for L = 40 and distances r < 20. The error 
results from a linear fit taking the statistical errors of 
the data into account. Taking the same raw data but ne- 
glecting their errors, yields At = 0.613±0.012. Note that 
the hypothesis (Qq) was formulated in "conformal units" 
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FIG. 9. Fluctuating values of the conductance g (or the 
transmission probability T) at the critical point of the net- 
work model. 1300 realizations at a fixed distance r = 10 are 
displayed. 
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FIG. 10. Scaling plot of the typical conductance for three 
different system sizes. The horizontal and vertical axes are 
logarithmic. 

(ro = 1) whereas the numerical data are represented in 
length units given by the lattice constant of the network. 
Therefore, tq is written explicitly in Eq. ( p2] ) whereas it 
is treated as unity in Eq. (|2l"|). 

Higher moments of InT can be obtained by using (T q ) 
as a generating function: 

(n = l + g(lnr) + ( g 2 /2)<(lnr) 2 ) + ... 

A straightforward calculation starting from ( po| ) yields 
for the log-variance (in conformal units) 

((In T) 2 ) - (In T) 2 = 16(/(0) - 4/'(0)) lnr 

n -2A M 8n y{\)d\ 
„ V ; (A 2 + l)cosh(^A/2) ' V ' 

The integral on the right-hand side decays algebraically 
to a constant whereas the first term grows logarithmi- 



cally. Eq. (|23|) says that the sum, E, of the log- variance 
and the integral is linear in lnr. To obtain the prefactor 
of lnr, we numerically calculated the integral for each 
value of r, took the log-variance from our data for sys- 
tems of size L = 60, and plotted the function £(r) versus 
X t \n(r/ro), see Fig. 11. It turned out that the integral 
was nonnegligible for the range of r values that were nu- 
merically accessible. The slope of the straight line in 
Fig. 11 is 2.08 ±0.11, i.e. the log-variance is about twice 
the log-average for large r (as is the case for in quasi-lD 
systems p2[). Given 8/(0) = A t , our numerical result 
implies that f'(0) is very close to zero: 

-4/'(0)/X t = 0.005 ± 0.008 , 

which leads us to conjecture /' = or, equivalently, a 
dependence of A\ which is exactly quadratic: 



A A = .,(0)(A 2 + 1)^(A 2 + 1) 



(24) 



Aside from being the simplest possible expression for A\ , 
this guess is in line with field theoretic expectations: in 
conformal field theories with a stress-energy tensor that 
is quadratic in the currents, the scaling dimensions are 
proportional to the quadratic Casimir invariant. The 
polynomial A 2 + 1 is, in fact, the quadratic Casimir of 
GL(2|2), evaluated on the continuous series of represen- 
tations A. In the language of multifractality (Sec. IV Ej ). 
the conjecture ( |24| ) means that the parabolic approxima- 
tion to the F(a) spectrum is exact. In the remainder 
of the present paper we shall assume the quadratic form 
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X t h(r/r ) 

FIG. 11. The log- variance of the conductance augmented 
by the integral on the right-hand side of (^) is denoted by 
E(r) and plotted versus -Xtln(r/ro). 



D. Multifractal spectrum 

We have seen that the dependence of (T q ) on the dis- 
tance between the two point contacts is governed in gen- 
eral by a continuous set of exponents 2Aa- This depen- 
dence simplifies, of course, in the asymptotic domain r — > 
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oo. For q > 1/2 the asymptotic behavior is controlled by 
the smallest exponent, 2Ao = X t /4 ~ 0.16. In the range 
— 1/2 < q < +1/2, the dominant contribution comes 
from the second term on the right-hand side of Eq. (20), 
with the exponent being 2A i ( 2g _i) = X t q(l — q) < X t /4. 
For q < —1/2 there appear additional contributions due 
to the poles of T(q-^±^) at A = ±i(2q+l), ±i(2q+3), 
etc. However, since A i(2g _i) < \( 2q +i) < \(2 q +3) < ■■■ 
for q < —1/2, these are negligible in the limit r — > oo. 
Thus we have 



conversely, a — X'(q(a)) = 0, F(a) = —aq(a) + X(q(a)). 
Using these relations we find 



(T q ) ' 
X{q) 



r -X{ q) j 

X t /4 for q > 1/2 , 
X t q(l - q) else . 



(25) 



Note that the spectrum of exponents X(q) is a nonde- 
creasing function of q. The spectrum shown in Fig. 12 is 
nonlinear, as is characteristic of a multifractal. 




FIG. 12. The multifractal spectrum X(q) governing the 
asymptotic power law scaling of the moments {T q ). 

That X(q) becomes constant for q > 1/2 can be traced 
to the boundedness of the variable InT < 0: for q > 1/2 
the asymptotics of the moment (T q ) is governed by the 
value of the distribution function at the upper bound, 
prob(T = l;r) ~ r~ x '/ 4 , which is independent of q, as 
is T 9 | T _ . In a thermodynamic interpretation of multi- 
fractal spectra, the nonanalyticity at q = 1/2 is called a 
"phase transition" in X(q) ff]. The failure of X(q) to 
become linear for q — > — oo results from the absence of a 
lower bound on InT. 

The multifractal nature of the X(q) spectrum can al- 
ternatively be described by a spectrum of f{a) type. If 
we make a multifractal ansatz for the probability density 
of T, 

prob(T = r~ a ) dT ~ r~ F(a) da , 
the moments of T for large r scale as 



-ia-F(a) da ^ r -X( q) 



Ha)= ( ^f (a>0). 



(26) 



The F(a) spectrum shown in Fig. 13 is defined only for 
positive values of a because a — X'(q) > 0. 



F(a) 




where X(q) and F(a) are Legendre transforms of each 
other, i.e. q + F'(a(qj) = 0, X(q) = qa(q) + F(a(q)), and 



FIG. 13. The multifractal spectrum F(a) describing the 
scaling of the distribution prob(T = r _a )dT. 

An obvious question arises: is it possible to relate the 
multifractal exponents X(q) to the r(q) spectrum of the 
critical eigenstates, or the scaling exponents of the LDoS? 
Our first observation is that there exists an obvious dif- 
ference between the two cases: while the spatial correla- 
tions of the LDoS continue to scale in a nontrivial manner 
with system size L when L is increased, the distribution 
of point-contact conductances (at fixed r) becomes inde- 
pendent of L |l5f| . Another way of saying this is that the 
critical conductance between interior point contacts has 
a trivial infinite volume limit, whereas the LDoS does 
not. 

To explain this distinction in physical terms, recall the 
iterative procedure by which the stationary limit ipoo was 
approached. An essential ingredient in that process was 
the draining action of the broken links at the two con- 
tacts. By the loss of probability through these outgoing 
channels, the system relaxes and settles down to a sta- 
tionary long-time limit. The key to understanding the 
size dependence of the conductance is to visualize the 
distribution of IV'ool 2 in space. In the regime of local- 
ized states this distribution is concentrated in a circular 
area of radius ~ £ centered around its source (i.e. the 
contact where the current is fed in). Outside this area 
the intensity I'i/'ool 2 falls off exponentially with distance. 
On approaching the critical point, £ diverges and the 
exponential decay turns into a power law. Thus the in- 
tensity becomes more spread out. However, in spite of 
this spreading out, the intensity remains (algebraically) 
localized near the source. The finite size of the network 
affects only the tails of the distribution l^ool 2 and, as a 
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result, the conductance converges to a well-defined limit 
as the system size goes to infinity. 

When the LDoS is discussed in the same language, a 
different picture emerges. Consider, for simplicity, the 
q th power of the density-density correlator 



1 



e-'En- 1 ^! 3 ' 



which, though not expressible as a correlation function 
of the LDoS, is related closely enough to allow a mean- 
ingful comparison. There again exists a dynamical inter- 
pretation. To compute the density-density correlator we 
feed in current through link c and, after relaxation to the 
steady state, measure the square of the amplitude at link 
d . What is different from the conductance measurement 
is that there are no drains at c and d (the network is now 
isolated). Instead, particles are absorbed at a constant 
rate 28 at any location in the network. As the system size 
is increased, the time spent in the network grows as the 
Heisenberg time <5 _1 ~ L d . (Previously, the dwell time 
in the 2D network was limited by the "strength" of the 
contacts and remained finite in the limit L — > oo.) Thus, 
there is a build up of particles in the network. When all 
states are localized, this causes the g-dependent correla- 
tor to diverge as d~ q ~ L qd . At the critical point, the 
divergence from build up is counteracted by the reduction 
in the amplitude at link j due to incipient derealization 
of the wavefunctions. For a general value of q the two 
competing effects do not cancel, so that a nontrivial scal- 
ing with L is expected to remain, in agreement with what 
we found in Sec. 0. 

Having said all this, we return to the original ques- 
tion: is it possible, after all, to relate the multifractal 
exponents X(q) to the r(q) spectrum? We wish to offer 
the following argument. Imagine placing tunnel barriers 
at the two contacts c, d. For zero tunneling probability, 
the network is closed and we can measure the density- 
density correlations. Now assume that the g-dependent 
density-density correlator has the same set of scaling ex- 
ponents as the LDoS correlator. Then, in the limit oj = 
(or L u — oo) we have 



(c 1(1 



• rto-'lc) 



2q 



As the tunneling rate is increased, a new time scale ap- 
pears: the time a particle injected at link c spends in the 
network (with the absorption rate 6 set to zero) before 
exiting through link d . The critical dynamics translates 
this dwell time into some characteristic length, L T . Al- 
though L T is an irrelevant length for an almost closed 
system with high tunnel barriers, it takes the regularizing 
role of the system size in the regime of open systems with 
L T < L. In the limit of vanishing tunnel barriers, the 
density-density correlator turns into the point-contact 
conductance. At the same time, the length L T must be- 
come proportional to the distance between the point con- 
tacts, for the simple reason that no other length scale re- 
mains available. This argument would say (T q ) ~ r~*( q ', 



and suggests z(q) as a candidate for X(q). We should be 
cautioned by the fact that the density-density correla- 
tor, unlike the conductance, does not respect any upper 
bound. This difference influences the tails of the distri- 
bution and changes the high moments, at least. In fact, 
we have shown X(q) to be constant for q > 1/2, whereas 
z(q) continues to increase. However, the tails of the dis- 
tribution should not affect the typical values, and there- 
fore one might expect X t = X'(0) = z'(0) = 2(a - 2). 
Recall that we found X t = 0.640 ± 0.009. Other groups 
quote values 2{a -2) = 0.56 ±0.04 @, 0.54 ±0.02 (e.g. 
||), and 0.58 ± 0.04 ^|. These values are not all mutu- 
ally consistent. We leave it as an open problem whether 
there is a flaw in the argument linking X(q) with z{q) or 
there is a real discrepancy. 



E. Quasi-ID limit 

We now endow the network model with a cylinder ge- 
ometry. This means that we consider an infinitely long 
strip of width W, with coordinates x £ K and y £ [0, W], 
and impose periodic boundary conditions in the trans- 
verse direction. As before, our interest is in the con- 
ductance between two point contacts, which are placed 
at positions (x,y) and (x',y'). The conductance in such 
a cylindrical setting can be related to the point-contact 
conductance in the infinite 2D plane with complex coor- 
dinate z by a conformal transformation 

2tt, 

z = exp W ' 

The conformal field theory rule for translating two-point 
functions from the plane to the cylinder reads EOj 



-2A* 



{ Vx {z)^{z')) = \z-z'\- 2 ^ 

W ■ , ( n , I • • A 

— Sinn I — [x — x + ly — ly ) 
7T V. W 

From this rule we have the relation 



irpq\ I cylinder _ , q . i2D 

\ >\(x,y);(x',y') \ ' lr=| ^ sinh(^- (x—x'+iy-iy')) | 

In particular, for the typical cylindrical conductance we 
obtain 



exp (In T) = 



w ■ 1 ( n , I ■ ls\ 

— smh \—{x-x +iy-iy)j 



(Recall that we are using length units so that (lnT) 2 D = 
for r = 1.) In the quasi-lD limit L = \x — x'\ 3> W, 
this result simplifies to 



exp (In T) 



(W/2w) Xt exp-irX t L/W . 



Note that from (|22|) the numerical value of irX t is 
7rA t = 2.01 ±0.03 . 



(27) 
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What is this result, a value of nX t close to 2, try- 
ing to tell us? Let us offer some speculation based on 
the assumption that the relation irX t = 2 holds exactly. 
In |fl7| it was shown that, if a naive continuum limit 
is assumed (i.e. possible rcnormalization effects due to 
short wave length modes are ignored) , the supersymmet- 
ric vertex model (|lj) for the critical network is equiv- 
alent to Pruisken's nonlinear a model |flf at couplings 
<3 X x — &xy — 1/2. The action functional of the latter 
model is 

S J dxdy (^(7 xx (^xx ~\~ £"yy) ~^ @ xy(£*xy £yx) 



where 



£ M „ = STr(l - ZZ)- 



ZZ)- x d v Z . 



For the sake of the argument, let us now assume Pruis- 
ken's model at a xx — o~ xy — 1/2 to be a fixed point of 
the renormalization group. Then, by raising the short 
distance cutoff from a = 1 to a = W we can reduce 
Pruisken's action S to a ID effective action 

W f ~ 1 -, 

Sid = — / dx STr(l - ZZ)- x d x Z(\ - ZZ)~ x d x Z . 

Alternatively, we could argue that in the quasi- ID limit 
the dependence of the field Z any can be neglected and, 
since a xx does not renormalize (by the fixed point hy- 
pothesis), the process of scaling out y simply produces a 
factor J dy = W. 

The one-dimensional theory with action functional 5id 
has been much studied, and its mean conductance is 
known [^2|-p[ to decay with length L of the conductor, 
which plays the role of distance between the contacts, as 



(T) 



ID 



exp-L/2W 



Moreover, from we know that the localization lengths 
for the mean and typical conductances of the one- 
dimensional nonlinear a model differ by a factor of 4, 
so that 



exp(lnT) 



ID 



exp-2L/W , 



which agrees, within the numerical errors, with what we 
found in ((2^). To turn the argument around, by assuming 
Pruisken's nonlinear a model at o~ xx = o- xy = 1/2 to be 
a fixed point, we would have predicted X t to be 



X t = {■KOxxY 



= 1/2 



0.637 . 



The above argument is not convincing, as it relies on 
the questionable assumption that the nonlinear a model 
is a fixed point theory. Conventional wisdom has it that 
critical two-dimensional nonlinear a models are unsta- 
ble with respect to quantum fluctuations and flow under 
renormalization to theories of the Wess-Zumino-Witten 
type. However, we can reformulate the argument and 



avoid any reference to Pruisken's theory. Let us assume 
that it is the network model itself (or, rather, a suitable 
continuum limit thereof) which is a fixed point of the 
rcnormalization group. Note that such an assumption is 
consistent with the fact that network model observables 
start scaling very rapidly when the observation scale is 
increased. (For example, in Fig. 10 there are no visible 
deviations from linearity as r approaches the short dis- 
tance cutoff a.) As before, we imagine raising the cutoff 
by using a sequence of RG transformations. By the fixed 
point hypothesis, we arrive for a = W at the ID net- 
work model (or, rather, some continuum version closely 
related to it), with the distance between contacts rescaled 
to L/W. Next we pass to the ID supersymmetric ver- 
tex model, and from there to the continuum action Sid- 
In contrast to earlier, the last step is benign, as the ID 
nonlinear a model is superrenormalizable (i.e. ultraviolet 
finite) and the RG trajectory can no longer depart from 
it. By this token, we again arrive at A t = 2/tt, this time 
without having passed through Pruisken's model. Thus, 
the proposed value for X t follows as a consequence of as- 
suming the network model (or a suitable continuum limit 
thereof) to be a RG fixed point. 

Note that if the above fixed point assumption is cor- 
rect and X t = 2/tt holds exactly, then we are led to the 
striking conclusion that Wess-Zumino-Witten models are 
ruled out as candidates for the fixed point theory. Indeed, 
the scaling dimensions for such models are given 46 by 
C\/(k + /i*), where C\ is the quadratic Casimir, and k 
and h* are integers. Such an expression for the scaling 
dimensions cannot produce the irrational number 2/tt. 



F. Reconstruction of the distribution 

We now return to the two-dimensional network and re- 
construct the entire distribution function for the point- 
contact conductance from the moments (T q ). By mak- 
ing a simple variable substitution (namely x = 2p+ 1) in 
Eq. (A4) of the appendix, the product of gamma func- 
tions in the formula for (T q ) can be represented as an 
integral over Legendre functions: 



\ni-\-%)\ 2 

T(qf 



(l + p)-«7> (i A-i)/a(2p + l)dp 



Next we define a probability density prob(p; r)dp for the 
variable p by 

/>oo 

prob(p; r) = / r- x <( A2+1 >/% A _ 1)/2 (2p + l) M (A)rfA . 
Jo 

Given 2A A = X t (X 2 + l)/4, comparison with Eqs. (|l|) 
and (|l^) yields 

/>oo 

(T q ) = / (l + p)- q W oh(p-r)dp . 
Jo 

Hence, on making the identification T = (l+p)^ 1 we con- 
clude that the probability density for T is prob(p; r)dp. 
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Although this is easily expressed in terms of T by using 
the inverse relation p = T^ 1 — 1, which has differential 
dp = — T~ 2 dT, we find it more convenient to work with 
the variable p instead of T. 

Because the Legendre functions "P(iA-i)/2(2p + 1) are 
oscillatory with respect to A (incidentally, they oscillate 
also w.r.t. p), the above formula for the probability den- 
sity is not well suited for numerical evaluation. Motivated 
by this, we switch to a different representation as follows. 
The Legendre functions satisfy the hypergeometric differ- 
ential equation 







d 



l( A + !) + TtAp + P«A-i)/a(2p + 1 



dp 



dp 



, 



and the integral of 7 ? (iA-i)/2 against the Plancherel mea- 
sure gives a Fourier representation of the ^-function: 



8(P)= / 7> ( a-i)/2(2p + l)MA)dA. 
Jo 

Both facts are standard results in harmonic analysis on 
the hyperbolic plane (or Lobachevsky plane) and are 
briefly reviewed in the appendix. Using them in the A- 
integral representation for prob(/3;r) we obtain 



^ r + x 4p p{p+i) ^)p Toh ^ r)=o > 

lim prob( j o; r) = 5{p) . (28) 

In r— >0 



Consider now the hyperbolic plane with the metric ten- 
sor in polar coordinates 9, (f> given by d6> 2 + sinh 2 (2<9)d(/> 2 . 
If we substitute l 



sinh 9, the differential operator 



d p p{p + l)d p turns into 



1 



— D— - - 

dp + 'dp ~ 4sinh(2<9) 89 



9 smh(20)^ 



which coincides with 1/4 times the radial part of the Lap- 
lace-Beltrami operator on the hyperbolic plane. There- 
fore, by viewing lnr as "time" and Xt/4 as a "diffusion 
constant" , we can interpret the initial value problem ( ps| ) 
as the heat (or diffusion) equation on that space. Solving 
the heat equation on the hyperbolic plane is a textbook 
example in Riemannian geometry p7j. For our purposes, 
a convenient expression for the solution is the following 
integral: 



prob(/o; r) 



27r -l/2 r -X t /4 



-t>/(X t lnr) tdt 



(Xtlnr)3/2 Kr v/sinh 2 W -/ 

arcsinn^/p y \ ) r 



which is easy to compute numerically. The result for the 
distribution function 



/ (In T;r) = prob(p; r/r ) 



dp 



d\nT 



(29) 



is plotted in Fig. 14 for the distance r = 15 between the 
contacts. The value X t — 0.64 is assumed. The error 



bars correspond to the mean deviation to be expected 
in histograms accumulated from 1760 independent mea- 
surements of data following the predicted distribution. It 
is seen that our analytical prediction agrees well with the 
numerical data points (accumulated from 1760 conduc- 
tances) represented by dots. 




In T 

FIG. 14. Normalized histograms of the critical log-con- 
ductance InT (dots) for distance r — 15. The numerical data 
are compared with plots of the distribution function /(In T; r) 
defined in (0). 



V. SUMMARY 

We have presented a numerical and analytical study of 
point-contact conductance distributions for the Chalker- 
Coddington network model. After reconsidering the mul- 
tifractal correlations of the local density of states, we 
first focussed on the distribution of point-contact conduc- 
tances in the quantum Hall plateau region, where strong 
localization of electrons occurs. As expected, the distri- 
bution is close to log-normal and is essentially parameter- 
ized by only the typical localization length. In particular, 
we found the log-variance to be proportional to the log- 
arithm of the typical conductance, with the constant of 
proportionality being —1.00 ± 0.05. 

We then turned to the plateau-to-plateau transition of 
the quantum Hall effect. Our analytical results are sum- 
marized as follows. By transforming the network model 
to a supersymmetric vertex model with GL(2|2) symme- 
try, we derived a formula, Eq. (|l7j), for the q th moment 
of the point-contact conductance at criticality. The gen- 
eral structure of the formula is completely determined by 
group symmetry. The unknowns are the scaling dimen- 
sions A a of certain local operators if\ q , which represent 
the point contacts in the formulation by the vertex model. 
We assumed these scaling dimensions to be proportional 
to the quadratic Casimir invariant of the symmetry alge- 
bra: Aa = X t (X 2 + l)/8. (This assumption is not essen- 
tial and can in principle be relaxed.) This choice leaves 
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X t as the only free parameter. Knowledge of all the mo- 
ments allowed us to reconstruct the entire distribution 
function. 

Salient predictions of our analysis are: i) The distri- 
bution of point-contact conductances becomes indepen- 
dent of the system size L in the thermodynamic limit 
L — ► oo. ii) At the critical point, the typical point- 
contact conductance of the infinite 2D network decays 
with the distance r between the two contacts as a pure 
power: exp(lnT) = (r/ro)~ Xt . iii) The log-variance 
equals —2 times the logarithm of the typical conduc- 
tance, iv) For large distances between the contacts, the 
q-moments of the conductance exhibit multifractal statis- 
tics: (T q ) ~ r ~ x{ -i\ where X(q) = -X t q(q - 1) for 
q < 1/2 and X(q) = X t /4 for q > 1/2. Thus there is 
a "phase transition" in the X(q) spectrum at q = 1/2. 

All these predictions are consistent with our numerical 
data, which were accumulated by a computing effort of 
about 2000 CPU hours on a Sun Sparc workstation. We 
found the distribution of point-contact conductances for 
r < L/2 to show no significant dependence on the sys- 
tem size, as expected. In a double logarithmic plot of 
the typical conductance versus r, the data points scatter 
around a linear curve with slope —X t — — 0.640 ± 0.009. 
The log-variance is linearly related to the logarithm of 
the typical conductance, with the constant of propor- 
tionality being —2.08 ±0.11. The phase transition in the 
X(q) spectrum is hard to see in our data, since the nu- 
merically accessible values of r are not large enough in 
order for the asymptotic behavior to dominate. However, 
the predicted distribution function for the point-contact 
conductances agrees well with our numerical data. 

On a speculative note we argued that, if the network 
model (or, rather, a suitable continuum limit thereof) is 
a fixed point of the renormalization group, then the scal- 
ing exponent for the typical point-contact conductance 
must have the value X t = 2/ir. This follows from confor- 
mal invariance linking 2D with quasi- ID, and from exact 
results available for the latter. While the fixed point 
assumption for the network model needs to be substan- 
tiated, it is remarkable that the predicted value lies very 
close to the numerical result. 



As a suggestion for further work, recall from Sec. [VE 
that conformal invariance at the critical point predicts 
the typical conductance between two point contacts at 
positions (0, 0) and (x, y) on a cylinder with circumfer- 
ence W to be 



cylinder 

(0,0) 



IT' 



sinh 



W 



[x + iy) 



Verification of this relation would provide a stringent test 
of the idea of a conformal fixed point theoy for the quan- 
tum Hall transition. We have not done the test, as our 
numerical calculations had already been long completed 
by the time we became aware of the exactness of the re- 
lation. We invite other groups to perform this stringent 
test and reduce the statistical error on X t . We feel cer- 
tain that the value of X t will be a benchmark for the 



analytical theory yet to be constructed, and is desirable 
to know with the same accuracy as the localization length 
exponent v. 
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APPENDIX A: 

Here we derive the result for the Clebsch-Gordan co- 
efficient and the Plancherel measure announced in equa- 
tions ([l9]) and (|l8|). For reasons that were explained in 
the text, this isn't an easy calculation. Fortunately, we 
can do it by using the following trick. 

We consider a very simple network, consisting of just 
two edges that interact along a chain of L vertices (see 
Fig. 15). 




FIG. 15. Graphical representation of a network, consisting 
of two edges that interact along a chain of vertices. 



The formalism developed in Sees. IV A and [VB ap- 



plies to this case just as well as to the two-dimensional 
network model. In particular, the q th moment of the 
conductance is given by a formula such as (^). A sim- 
plifying feature is that the product of 1Z matrices now 
organizes into a convolution product of transfer matri- 
ces, T L . Denoting the eigenvalues of the transfer matrix 
T by t(X) we get 



(T q ) 



1D:L — 



(Vq, V*0\Xq) \ t(\) L /i(\)d\ , (Al) 



by a similar reasoning as in the body of the paper. Our 
strategy will now be to exploit the simplicity of this ID 
model and compute the moments (T 9 )id ; l from a quite 
differe nt a pproach. By comparing the result to the for- 
mula ( |Al[ ), we will ultimately be able to read off the de- 
sired expressions for the Clebsch-Gordan coefficient and 
the Plancherel measure. 

For technical convenience, we shall consider the two- 
edge network model in the limit of weak backscattering 
at the nodes. An attractive feature of this limit is that 
the task of computing (T 9 )id ; l can be reformulated as a 
partial differential equation (of the Fokker-Planck type) 
which is readily solved. Before writing down that equa- 
tion, it is helpful to make two adjustments. Rather than 
computing directly the moments (T q ), we will study the 
entire distribution function of T. Also, we switch from T 
to the variable p = (1 — T)/T ("Landauer's resistance"). 
Now, by an elementary calculation (for a review, see e.g. 
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|]48[|), the probability density prob(p; L)dp of Landauer's 
resistance satisfies the differential equation 

e— W ob{p- L) = —p(p + 1)— prob(p; L) (A2) 

where I is the elastic mean free path. We are going to 
solve this equation by harmonic analysis, i.e. by diago- 
nalization of the differential operator d p p(p+ l)d p . [This 
operator has a geometric meaning as the radial part of 
the Laplacian on a noncompact Riemannian symmetric 
space SU(1, 1)/U(1).] Introducing the Legendre function 
V v {x) through its integral representation, 

V v (x) = J J ^ (x + y/x 2 -l cos 0) " , (A3) 
one easily verifies 

§~ p P(P + 1)^ + \(>? + !)) n^i)/ 2 (2p + 1) = . 

This relation suggests a solution of the differential equa- 
tion (A2) of the form 



POO 

probfc; L) = / e-6< Aa+1 >P (4A _ 1)/2 (2p + l)dm(X) . 
Jo 

The spectral measure (or Plancherel measure) dm(X) is 
determined by the asymptotic behavior of the Legendre 
functions 'P(iA-i)/2(2p + 1) for p — > 00, as follows. By 
using the su bsti tution u = tan(</>/2) in the integral rep- 
resentation (A3), one finds 

7> (i A-i)/ 2 (cosh20) 9 -=3? c- e (c(\y xe + c(-A)e-* Ae ) 
where the c-function is given by 

1 r(|»A) 



c(A) 



/5Fr(i(<A + l)) 



From this asymptotic limit, we infer the orthogonality 
relations 



P(u-i)/2(2p+l)7V_i) /2 (2p+l)dp = N X S(X-X') , 



with the normalization factor being N\ — 7r|c(A)| 2 . In 
conjunction with the initial condition 

prob(p; L) = S(p) , 

stating that transmission through a short chain is ideal, 
the presence of the normalization factor determines the 
spectral measure to be 

7 /■> n dA X . /ttA\ 

dm(A) = ^ATF = 2 tan H^' 



To recover the moments (T 9 ) from the distribution 
prob(/?; L)dp, we need the integral 



(l+p)- q V (lX -i)/2(2p+l)dp, 



which converges for q > 1/2. We claim that this integral 
has the value 



h(q) = 



r( g -i + f)r( g -i-f 
r( q y 



To prove this statement, we proceed as follows. In the 
first step, we set x = 2p + 1 and write 



(x + l)-*P liX _ 1)/2 (x)dx 



= T(q)- 1 £ (J™ e - ax P (iX _ 1)/2 (x)dx\ a^e^da 

which is a valid equality for q > 0. According to Ref. |4! 
(p. 323, no. 11) the integral in parentheses equals 



c V(i\-i)/2(x)dx 



Ki\/2(a) 



This leads to an integral over the auxiliary variable a, 
which converges for q > 1/2 and the value of which we 
take from Ref. M [p. 716, no. 6.628(7)]: 



i"- 3 l 2 e- a K iX/2 {a)da 



r( g -l-^)r (g -l + ^)ii m ^-x)/ a (^) > 

KH 2 2' W 2 2 a^o (sinha)"?- 1 

where V£(x) is the associated Legendre function. From 
Ref. |3 (p. 332, no. 8.1.2) this function has the follow- 
ing small-a limit: 

^ 1)/2 (cosha) ^° r( g )- 1 (a/2)'- 1 . 

Combination of all these results yields 

h (g) = 2"- 1 f (1 + z)- 9 7V-i)/2 



(A4) 



r( g -i-f)r( g -i + f) 
r(g)a 



which proves the claim. 

By substituting T = (1 + p) _1 and inserting for the 
probability density prob(p; L)dp the spectral resolution 
given earlier, we finally arrive at 



1D-.L — 



r^ 2+1 ^/ u h(q)dm(X) . 



Comparison with (Al) identifies the eigenvalue of the 
transfer matrix T as t(A) = q~( x2+1 )/ ai , and yields 
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I (Vq, V*0\Xq)\ fi(X)d\ = h{q)dm(X) 



(A5) 



Although this result gives an answer for the product of the 
squared Clebsch-Gordan coefficient with the Plancherel 
measure, it does not allow to make the separate identifi- 
cations proposed in Eqs. ( |l9| ) and jig)). For that, more 
detailed considerations are necessary. For brevity, we re- 
frain from elaborating on these since, actually, all that is 
needed for the main text is the formula (pV5|). 
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